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This paper studies the impact of incorporating wind power generation WPG 
on the power system on prsence of voltage source converter based high 
voltage DC (VSC-HVDC). A new meta-heuristic optimization technique are 
use for solving of the optimal power flow (OPF) problem, this technique 
optimization namely Ant Lion Optimizer (ALO). The optimization method is 
the Ant Lion Optimizer (ALO) method for resolve the optimal power flow 
(OPF) with incorporating of wind power generation on prsence of VSC- 
HVDC. And we used weibull distribution model of the wind farm. The ALO- 
OPF method has been examined and tested on standard test systems IEEE 30 
bus with objective functions is minimization of cost total of production TPC 
are contain the sum of thermal and wind generation cost. 
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1. INTRODUCTION 

In recent years, consumption of electric power has increased significantly, so it is necessary to find a 
solution to this increasing demand. Therefore, the incorporating of renewable energy as wind energy into the 
energy system is one of the best solutions to solve this problem. This source is considered one of the 
desirable renewable energy sources because it is not polluted and not expensive. Integrating wind farms in 
presence of High Voltage Direct Current (HVDC) to power grid will inevitably present a big impact on the 
electrical power system safety and economical operation of the system and such impact has become a lot 
research topic as wind farm power penetrations increase in power system in recent years [1]. The VSC- 
HVDC transmission is gaining a lot of interest in the power industry, and are much more preferable than AC 
ones for high voltage levels because they are more economic and more reliable, transmission capacity is 
greater in HVDC lines due to the non-existence of the skin effect, and also the efficiency and controllability 
of DC converters are higher [2]. In this scenario, it necessary to solve the optimal power flow (OPF) problem 
by using Meta heuristic algorithms when the wind power generation (WPG) are introduced in electrical 
system [3] in presence of VSC-HVDC, where the objective is to find the best operation of electrical power 
system which minimizes the total production cost and reduction of power losses for a best balance of 
electrical power system. 

Numerous Meta heuristic algorithms are used for solving OPF problem. Several methods are 
implemented on optimal power flow problem such as, differential algorithm (DE) [4,5], Grey wolf optimizer 
(GWO) [5], moth-flame optimizer (MFO) [6], moth swarm algorithm (MSA) [7],artificial bee colony 
algorithm (ABC) [8] and particle swarm optimization (PSO) [9], The Ant Lion Optimization (ALO) 
technique is one of the modern meta-heuristic algorithms for solving the optimization problem, more 
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research applied this algorithm for solving the optimization problem such as in[10], the authors applied ALO 
algorithm for solving optimal power flow problem with objective function is minimization of total fuel cost, 
total gas emission and total active losses. Also in [11], the authors applied ALO for objective functions are 
minimization of Fuel Cost, Voltage Deviation, and improvement of Voltage Stability. 

Nowadays, the optimal power flow for the system that includes wind farm power is the subject of 
ongoing research model. In this paper optimal power flow is carried out for combined wind power farm to 
electrical power system with and without in presence of VSC-HVDC by using the ant lion optimization 
(ALO) technique and using weibull distribution model of the wind farm, is best suited probability density 
function for modelling the wind farm output [12]. Panda and Tripathy in [13] applied modified bacteria 
foraging algorithm (MBFA) for solving optimal power flow of wind integrated power system with an 
objective functions is minimization of total cost production (thermal and wind generation) and in[14] applied 
MBFA for solving security-constrained optimal power flow solution of wind-thermal generation system. In 
[15], the authors applied three powerful meta-heuristic algorithms-namely, the cuckoo search algorithm 
(CSA), firefly algorithm (FFA), and flower pollination algorithm (FPA) to optimize power flow in Algeria's 
Adrar power with Multi-objective function includes: the fuel cost, the NOx emissions, and the imbalance cost 
of the WPGs. Kiran et al [16] applied Multi-objective hybrid PSO-APO algorithm to solve the multi- 
objective security constrained optimal power flow (MO-SCOPF) problem. 

The objective of this paper is to optimize an objective function to solve constraint based quadratic 
cost function. Ant Lion Optimization is employed for solving an optimal power flow problem includes wind 
power in presence of VSC-HVDC. The proposed algorithm is applied to minimization of cost total of 
production contain the thermal and wind generation cost. The proposed algorithm is tested on IEEE 30-bus, 
and results of simulation are compared with two methods WOA and GWO methods. 


2. MODEL OF VSC-HVDC 

This section introduces the modeling of the AC/HVDC system. The modeling is based in the 
converter station and DC grids system models as shown in Figure 1. The model a converter station consists 
of a controllable voltage source U, = U,Z6,behind the phase reactor, represented as complex impedance 
Z, = R, + jX¢. the low pass filter from Figure 2 is represented as a susceptance Brat system frequencies. A 
transformer connects the filter bus to the AC grid and is represented by its complex impedance Z;¢ = Ry + 


jX (171. 
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Figure 1. VSC MTDC converter station [17] Figure 2. Equivalent single phase power flow 
model 


The equations for active and reactive power at the grid side in terms of the complex voltages are: 

Ps = —U3G;¢ + UU; [Grp cos(d, — 5¢) + Bey sin(5, — 6)| (1) 

Q, = U2Bir + B.By [Ger sin(S, — 5) — Br cos(6, — 5¢)] (2) 
The equations at the converter side are: 

P, = U2G, — U;U,[G, cos(5; — 6.) — B, sin(5; — 6,)] (3) 

Q, = —U2B, + U;U,[G, sin(6; — 5.) + B, cos(6; — 5,)] (4) 


Indonesian J Elec Eng & Comp Sci, Vol. 12, No. 2, November 2018 : 625 — 633 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 0 627 





The expressions for the filter side complex power flowing through the transformer are written as 
Pop = U? Gee — UUs [Grp cos(5, — 5) — Bry sin(5, — 5¢)| (5) 


Qs = —Uf Ber + UpUs|Gep sin(5s — 57) + Ber cos(5s — 6/)] xe 
And those flowing through the phase reactor side are 
Pog = —UPG, + UrU.[G, cos(5; — 5,) + B, sin(5; — 6.)] “2 
Qcr = UB, + UrUc|G, sin(5; — 6.) — B, cos(5- — 6.) © 
The DC system modeling can be represented by a resistive network with current injections and DC 


voltages at the different nodes[18]. The current injected at a DC node i and the power flow equation of DC 
line can be written as: 


lac; ~ Yacij . (Vac = Vac,) (9) 


Pac; = PU ac; Yacij: (Vee on Uae;) (10) 


where Yaci; equal to 1 / Racij and p =lfor a monopolair system or p =2 for a bipolair system. 


3. FORMULATION PROBLEM 
3.1. Problem 

The objective of the OPF problem is to minimize the total cost of all available generators plants 
subjected to the operating constraints (equality and inequality constraints) in the power system. In this paper, 
the objective function is the sum of the operating costs of each available conventional generator and the wind 
farms can be expressed as follow: 


F= yee C.(Pei) + L¥pws (PW as _ Pw;) + Y Im, ws (Pw; — Pwjav) (11) 


In this expression. The first term denote the cost of thermal power generation, second term is the 
cost due to under estimation of available wind power and third term is the cost due to over estimation of 
available wind power, Details of the first, second and third terms are explained as: 


3.2. Fuel cost of the conventional generator 
Widely acclaimed quadratic cost function for the thermal generators are mathematically represented 
by (12) as follows: 


C,(Pgi) = Ay + ByPai + CPG; (12) 
The cost of wind power generation (WPG) is calculated considering the over-generation and under- 
generation of the available wind power, based on the Weibull distribution function WDF of the wind speed 


[19], represented by the second and third terms of Equation (11). 


3.3. Cost function due to the over-generation 
The cost function due to over generation of available wind power is expressed as follow [19]: 


CowiCWiav _ W;) = ky i (Wie _ Wi) -_ Ry In (W a Wi fwW) (13) 


Where k,; is the penalty cost coefficient for over-generation of if? wind farm, fy(W) is the probability 
density function (PDF) of wind power output. 


3.4. Cost function due to the under-generation 
The cost function of i‘” wind farm for calling the reserves to cover i*’ wind farm due to under- 
generation is expressed as follows [20, 21]: 
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Where k,,; is the reserve cost coefficient for under-generation of i” wind farm. 
The probability density function (PDF) of the Wiebull distribution is given by: 


1 (2 fis 


Where v is the wind speed, k and c are the shape factor and scale parameter (m/s). 
For the Weibull distribution function, the discrete portions of wind energy conversation system (WECS) 
power output random variable as following [19]: 


Pratea{W = 0} = 1 —exp(—(“2#)') + exp (- ((#=4)')) (16) 


fv) =*(2)" 


Cc 


Cc 


Pratea\W = W,atea} = exp (- (“eee ") + exp (- ((2=24)')) (17) 


Cc 


While the continuous portion of WECS power output random variable in becomes: 


klvcut—in (A+e*1)¥cut-in am (1+p*1)Vcut—in ‘; 
fw(W) = etch (eee Peue— in) exp (- (eer eur=in ) (18) 


Cc Cc 


4. THE ANT LION OPTIMIZATION ALGPRITHM (ALO) 

Ant lion optimization algorithm ALO is a new nature inspired proposed by Seyedali Mirjalili in 
2015 [22], for to solve any optimization problems of constrained engineering. The ALO algorithm modelled 
based on hunting mechanism of antlions in nature. The benefits of ALO algorithm, it can find optimal 
solution for minimizing the objective function. In ALO algorithm five main steps of hunting prey as the 
random walk of ants, trapping in antlions traps, building traps, sliding ants against toward antlion, Catching 
preys and rebuilding traps are implemented. 

The ALO algorithm mimics the hunting behavior of ant lions, i.e., the interaction between predator 
(ant lions) and prey (ant). 

In this algorithm, the expression mathematically of is the random walks of ants for detects the 
location of food is describes in equation (19). 


X(t) = [0, cumsu(2r(t,) — 1), cumsu(2r(t2) — 1),...,cumsu(2r(t,,) — 1)] (19) 


where X(t) is the random walks of ants, n is the max-iterations, t is the step of random walk, and r(t) is a 
function defined as follows: 


_ (lifrand > 0.5 
r= {9 if rand < 0.5 (20) 
where, rand is a randomly generated number uniformly distributed in the range of [0,1]. 

The Details of different steps of ALO algorithm describe the relation between predators (ant lions) 
and prey (ant) are explained as. 


4.1. Random walk of ants 

In every step of optimization, ants update their positions with random walk and Random walks are 
all based on the Equation (21), all the positions of ants are inside the boundary of the search space, since 
every search space has a boundary, the random walks of ants is normalized by using the following equation: 


si-a)eoted 
xf = GedPined 5 of (21) 
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where the a,; b; are respectively the minimum and maximum of random walk corresponding of i” variable. 
cf; df: are respectively indicated the minimum and maximum of ith variables at tth iteration. 


4.2. Trapping in antlions traps 
The influence of antlions traps on random walks of ants is expressed as follow: 


cj = Antlion; + c (22) 
dj = Antlion; + d‘ (23) 


4.3. Building traps 
During optimization, The ALO algorithm is required to use a roulette wheel operator for selecting 
antlions based of their fitness for give a high chance for catching ants. 


4.4. Sliding ants against toward antlion 

With the mechanisms proposed so far, antlions are capable to build traps proportional to their fitness 
and ants are required to move randomly near of the center of pit. However, once antlions realize that an ant is 
in the trap, they shoot sands outwards the center of the pit. This mechanism mathematically is proposed 
modeled as follow: 


cha (24) 


dé =— (25) 


4.5. Catching preys and rebuilding the traps 

The final step of hunting is when an ant reaches in the bottom of the antlion pit and is caught in the 
antlion’s jaw. After this step, the predator pulls the ant inside the sand and consumes its body. Then the 
antlion is required to update its position to the latest position of the hunted ant to enhance its chance of 
catching new prey. Catching prey and re-building the pits calculated as given by the following equation: 


Antlion; = Antj,if f (Ant{) > f(Antlion;) (26) 
where Antlion{, Ant; is jth the position of the selected antlion and ant at iteration t. 


4.6. Elitism 

Elitism is one of the most important characteristic of all evolutionary algorithms. Elitism of ant lion 
is calculated by using the roulette wheel. In this study, the best solution(s) obtained by the ALO algorithm at 
any iteration is saved as elite. The elitism mechanism mathematically described as follows equation: 


t t 
Ant! = “Ate (27) 


Where R§ A is the random walk around the antlion and is selected by using the roulette wheel at tth iteration, 
RE Eis the random walk around the elite at tth iteration 


5. RESULTS AND ANALYSIS 

In this study for the purpose of verifying the performance and efficiency of the ALO algorithm 
proposed on the power system with and without incorporating wind power generation and presence of VSC- 
HVDC, tests are carried out on IEEE 30 bus system.. All the simulations are carried out by using MATLAB 
2009b, and computed Core (TM i5; CPU 1.80 GHz-8Go RAM). For establishing the superiority of the 
proposed ALO, 30 independent trial runs are performed for all the test cases with a comparative study 
reported in the following section, the population size (NP = 40) and number of iterations maximal (iter Max 
= 500) are taken same for three optimization methods ALO, GWO and WOA. 


5.1. Results of IEEE 30-bus test power system without WPG and VSC-HVDC 
To verify the proposed approach and for comparison purposes, we perform simulations on the IEEE 
30-bus system, involves 6 generators, 41 lines, 4 transformers. About this case study, 9 reactive 
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compensators are installed at buses. The load demand is 283.4 MW and 73.5 MVar n 100 MVA base. The 
results presented in Table | including the generation cost and real power losses of the proposed method and 
with other existing optimization methods. 


Table 1. Comparison of Fuel Cost Obtained By Proposed Algorithm and Other Algorithms. 
Variables min max ALO GWO WOA  DE(4) GWO(5) MFO(6) MSA(7) 


Pigss (MW) 8.6026 8.8347 8.7405 8.6150 9.30 9.3000 9.0345 
Cost ($/hr) 799.0133 799.2954 799.8095 799.2891 801.41 799.072 800.5099 














In order to verify and test the robustness and performs of the proposed algorithm for optimal 
solution of minimizing the objective function, we made a comparison between of the proposed method with 
other existing optimization methods presented in the literature. Here, we considered fuel cost of thermal 
generators as an objective function. The comparison of minimum fuel cost and minimum power losses 
obtained by the proposed method (ALO) with other existing methods as GWO, WOA, DE, GWO, MFO, and 
MSA has been presented in Table 1. From the results, the minimum fuel cost and minimum power active 
losses obtained by the proposed ALO method are 799.0133 $/h and 8.6026 MW respectively, which are 
better than all other algorithms as well as previously reported methods. The proposed ALO optimization 
method superior and robust compared to GWO and WOA algorithms in terms of getting near global optimum 
point and faster convergence for solving the OPF problem. The convergence characteristics of the proposed 
method and, GWO and WOA methods are shown in Figure 3. 
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Figure 3. IEEE 30-bus system: Convergence plot for ALO, GWO and WOA methods 


5.2. Results of IEEE 30-bus with WPG and VSC-HVDC 

In this cases, we have study the influence of wind power generation WPG on power system in 
presence of VSC-HVDC device as minimization of total production cost (TPC) and reduction of power 
losses, The proposed ALO algorithm is tested by using Weibull distribution model of the wind farm with 
objective function is minimization of thermal generation cost (TGC) and wind generation cost (WGC), and 
its performance is compared with the standard GWO and WOA algorithms. To verify the performance and 
effectiveness of the proposed algorithm of solving OPF problem including wind power farm with and 
without presence of VSC-HVDC, two different case studies are considered, the case | in this studies is the 
base case (including wind power farm without VSC-HVDC), the second case is including wind power farm 
with presence in VSC-HVDC. Two wind farm (WF) consists of 43 units of wind turbine generation (WTG) 
are connected to the system at bus no 10 and bus 24, a nominal power rating of each WTG is 2 MW with a 
total capacity are 85.02 MW, the wind power generation of WF1 at bus No. 10 having a total capacity of 
56.68 MW and the wind power generation of WF2 consists a total capacity of 28.34 MW, This two values of 
wind power was taken from [1] and are considered as the best penetration of wind power. 

The parameters of wind turbine are the factors of wind speed distribution are: shape factor, k = 1.12 
and scale factor, c = 5.92, and the characteristics wind turbine on this studies are: P,aqteq = 2000 KW, 
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Vratea = 10.28 m/s, Veut-in = 4™M/S, Veut-o¢¢ = 25 m/s. The penalty factor is k, ; = 0.03 $/MWh and 
the reserve factor is k,; = 0.03 $/MW. 


Case 1: ALO Solution without VSC-HVDC 

The first case studies, the proposed ALO algorithm are tested for case of increased load at 45 % 
from total load for see the effect of wind energy on the system in case of contingency, the total load was 
410.93 MW. The results obtained by the proposed method and with those found by GWO and WOA are 
presented in the Table 2. 


Table 2. Comparison of fuel cost obtained by proposed algorithm and other reported algorithms case of 








contingency 
Variables Without WPG ALO GWO WOA 
Pross (MW) 13.7979 11.4941 11.7777 11.7524 
Thermal generation cost ($/hr) 1339.6 970.8575 971.1843 977.8275 
Wind generation cost ($/hr) /II 19.9495 19.9529 19.9515 
Total production cost($/hr) 1339.6 990.8070 991.1372 997.7790 





The Table 2 shows, real power losses and total production cost (TPC) by the proposed method and 
two optimization methods in literature are GWO and WOA. The active power loss was reduced from the 
value of 13.7979 MW (without WPG) to 11.5420 MW in case of integration of WPG. Therefore, it is clear to 
see effect of wind power generation (WPG) on the power system in case of contingency. The results obtained 
by ALO algorithm shown better results for the real power losses and the total cost of production in the values 
11.5420 MW and 990.0337 $/hr respectively compared results by GWO and WOA. 


Case 2: ALO Solution for over-load with WPG and VSC-HVDC 

In this part, we have use VSC-HVDC device for a best soultion in problem of over-load, transprted 
maximum energy and reduction of active power losses. In this case study, a VSC-HVDC is integrated in the 
IEEE 30 bus system between bus 2 and 5. The DC system voltage is +/- 132 kV. The converter station VSC1 
(slack bus of the DC network) is connected to bus 2 and VSC2 is connected to bus 5 of the AC system. The 
parameters and at the converter stations are given in Table 3 and the results obtained by the proposed method 
are presented in the Table 4. 


Table 3. VSC Converter Data 


Rating & Converter loss data 








Converter parameters 








No. 1,2 
Xt (p-U) 0.1121 t+Vac 132 
R,, (p.u) 0.0015 P..(MW) 90 
By (p.u) 0.0887 a 1.103 
X, (p.u) 0.1642 b 0.887 
R, (p.u) 0.0001 Cree 2.885 
Raco-s(P-U) 0.0472 ae 4371 





Table 4. Comparison of results obtained by proposed algorithm with and without VSC-HVDC 








Variables Without WPG and HVDC With WPG With WF and HVDC 
Pross (MW) 13.7979 11.4941 11.0183 
Thermal generation cost ($/hr) 1339.6 971.1843 970.1096 
Wind generation cost ($/hr) 19.9495 19.9529 19.9529 
Total production cost$/hr) 1339.6 991.1372 990.0625 





from Table 4, the presence in HVDC line between buses 2 and 5 reduce the total production cost, 
also the active power losses was reduced from the value of 11.4941MW (without HVDC line) to 11.0183 
[MW] in presence of HVDC line. Also from Figure 4, after the increase in the load at 45% (case without 
presence in HVDC line ), the power transmitted in line 5 between bus 2 and 5 is 80.00 MW (this value is the 
maximum limit 80 MW) and the power transmitted in line 1 between bus 1 and 2 is 127.7722 MW (this value 
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is near the maximum limit 130 MW), after integration of HVDC in line 5 (line 5: is optimal placement of 
HVDC device in IEEE 30 bus), the power transmitted in line 1 is decreased to equal 125.33 MW, the power 
transmited in HVDC line is 90 MW, and see that also the effect of VSC-HVDC to find a best solution for the 
problem of over-load are shown in Figure 4 by using the proposed method. 
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Figure 4. The effect of WPG on the Active powers transit with and without HVDC Line 


6. CONCLUSION 

In this paper, a recent metaheuristic technique namely Ant Lion Optimizer (ALO) was successfully 
applied to solve the optimal power flow (OPF) problem for case of integration of power wind farm on the 
power system with and without in presence of VSC-HVDC. We studied an optimal power flow for 
minimization of total production cost. The proposed ALO method has been examined and tested on the IEEE 
30-bus. The obtained results of ALO method were compared with those obtained from two methods 
developed by us, namely, Grey wolf Optimization and the Whole Optimizer Algorithm. The Ant Lion 
Optimizer (ALO) method has been a best result. The Simulation shows the best results for minimize the cost 
of generation, the active power losses and reduction of over-load in the transmission line when the WPG 
includes in the power system in presence of VSC-HVDC. 
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